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Abstract We present a cluster analysis of the bright main-sequence and faint pre- 
main-sequence stellar populations of a field ~ 90 x 90 pc centered on the Hll region 
NGC 346/N66 in the Small Magellanic Cloud, from imaging with HST/ACS. We 
extend our earlier analysis on the stellar cluster population in the region to char- 
^ c*j acterize the structuring behavior of young stars in the region as a whole with the 

Oh use of stellar density maps interpreted through techniques designed for the study of 

the ISM structuring. In particular, we demonstrate with Cartwrigth & Whitworth's 
parameter, dendrograms, and the A -variance wavelet transform technique that 
the young stellar populations in the region NGC 346/N66 are highly hierarchically 
clustered, in agreement with other regions in the Magellanic Clouds observed with 
HST. The origin of this hierarchy is currently under investigation. 



1 Method: The identification of stellar clusters 

For the investigation of the clustering behavior of stars it is necessary to thoroughly 
characterize distinct concentrations of stars, which can only be achieved by the ac- 
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Fig. 1 Schematic of the star-count process, (a) The chart of an observed stellar sample, (b) The 
corresponding stellar density map, after counting stars in quadrilateral grid of elements (pixels) 
of size 1.8" each, and after filtering the map with a Gaussian of FWHM ~ 2.8 px (~ 5"). (c) The 
corresponding isodensity contour map. Isopleths at levels ^ 3<7 are indicated with white lines. 



curate identification of individual stellar clusters. Considering the importance of this 
process, different identification methods were developed, which can be classified in 
two families. The first, represented by friend of friend algorithms and cluster analy- 
sis techniques, e.g., [2], are designed for limited samples of observed stars, and thus 
are based on linking individual stars into coherent stellar groups. These methods 
are recently superseded by minimum spanning trees, e.g., [1]. The second family of 
identification codes, represented by nearest-neighbors and star-counts, make use of 
surface stellar density maps constructed from rich observed stellar samples. Distinct 
stellar systems are identified as statistically significant over-densities in respect to 
the average stellar density in the observed regions, e.g., [5]. Tests on artificial clus- 
ters of various density gradients and shapes showed that the latter (density) tech- 
niques are more robust in detecting real stellar concentrations, provided that rich 
stellar samples are available [10]. A schematic representation of stellar density maps 
constructed with star-counts is shown in Fig. 1. 



2 Data: Stellar clustering in the region NGC 346/N66 

One of the most prominent bright stellar systems in the Small Magellanic Cloud 
(SMC) is the stellar association NGC 346, related to the Hll region LHA 115-N66 
[6], the brightest in this galaxy. This system appears in partially-resolved observa- 
tions form the ground as a single stellar concentration, but recent imaging with the 
Advanced Camera for Surveys onboard the Hubble Space Telescope (HST) allowed 
the detection of smaller sub-clusters within the boundaries of the Hll nebula. The 
images were collected within the HST GO Program 10248 and were retrieved from 
the HST Data Archive. Their photometry demonstrated that the faint young stellar 
populations in the region are still in their pre-main-sequence (PMS) phase, and re- 
vealed a plethora of sub-solar PMS stars [4]. Our nearest- neighbor cluster analysis 
of the observed young stellar populations, i.e., the bright main-sequence (down to 



Hierarchically Clustered Star Formation in the Magellanic Clouds 



3 



Fig. 2 Isodensity contour 
map from star-counts of the 
young bright main-sequence 
and faint PMS populations 
identified with HST/ACS in 
the region of NGC 346 in the 
SMC. Lines represent iso- 
pleths of significance ^ lo~. 
Apart from the dominating 
central large stellar aggregate, 
there are peripheral young 
sub-clusters, revealed as sta- 
tistically important stellar 
concentrations. The central 
aggregate, denoted by the 1 a 
isopleth, encompass various 
distinct sub-groups, which ap- 
pear at higher density thresh- 
olds. NGC 346 itself appears 
at ^ 3(7 significance. 



m 555 ~ 21) and the faint PMS stars, revealed a significant number of smaller, previ- 
ously unresolved, young stellar sub-clusters [9]. This clustering behavior of young 
stars in NGC 346 is further demonstrated here by the stellar density contour map of 
Fig 2, constructed with star-counts. 



3 Results: Hierarchical clustering of young stars 

The map of Fig. 2 shows significant sub-structure, in particular within the 1 a bound- 
aries of the central dominant stellar aggregate. This structuring behavior indicates 
hierarchy. The minimum spanning tree (MST) of the young stars in the whole region 
allows to determine the statistical parameter, introduced by [3] . This parameter is 
a measure of the fractal dimension D of a stellar group, permitting to distinguish be- 
tween centrally concentrated clusters and hierarchical clusters with fractal substruc- 
ture. The application of the MST to our data shows that the region NGC 346/N66 is 
highly hierarchical with a J2 that corresponds to a fractal dimension D ~ 2.5. 

Constructing surface stellar density maps allows us to further characterize the 
clustering behavior of stars with the application of tools, which are originally de- 
signed for the study of the structuring of the interstellar medium (ISM), as observed 
at far-infrared or longer wavelengths. The so-called dendrograms are used for the 
visualization of hierarchy through structural trees [8]. The dendrogram of the stellar 
density map of NGC 346 demonstrates that the observed hierarchy is mostly due to 
the substructure in the dominant stellar aggregate. The A -variance analysis [11, 7] 
is a robust structure analysis method that measures the amount of structure on a 
given scale I. In principle the A -variance is directly related to the power spectrum of 
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Fig. 3 The A -variance spec- 
tram of the surface stellar 
density map of the entire 
region of NGC 346/N66. 
This analysis shows that the 
young stellar populations in 
this region are hierarchically 
structured up to length-scales 
of ~40". The spectral index 
j8 is determined from the fit of 
the spectrum for data between 
lags 4" and 13"(indicated by 
the gray shaded area). The 
dashed line provides the used 
virtual beamsize (5"). 

the map, and thus for a power law spectrum of index — J3, A -variance also follows 
a power law, a\ « l a , with a = /3 — 2. The application of the A -variance analysis 
on the surface stellar density map of NGC 346 verifies that indeed the clustering of 
the young stars in the region is self-similar (Fig. 3), with a spectral index j3 ~ 2.8, 
corresponding to a fractal dimension D = 2.6 of the corresponding fractional Brow- 
nian motion structure [11], similar to that previously derived for Galactic molecular 
clouds. Self-similarity appears to brake, i.e., we find different hierarchical proper- 
ties for the short-range scaling and the behavior at the overall scale of the region, 
at length-scales I > 25 px, corresponding to physical scales of ~ 40"(~ 1 1 pc at the 
distance of the SMC). 
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